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ABSTRACT 

We investigate the characteristic of microlensing signals of Earth-like moons orbiting ice-giant planets. From 
this, we find that non-negligible satellite signals occur when the planet-moon separation is similar to or greater 
than the Einstein radius of the planet. We find that the satellite signal does not diminish with the increase of 
the planet-moon separation beyond the Einstein radius of the planet unlike the planetary signal which vanishes 
when the planet is located well beyond the Einstein radius of the star. We also find that the satellite signal 
tends to have the same sign as that of the planetary signal. These tendencies are caused by the lensing effect 
of the star on the moon in addition to the effect of the planet. We determine the range of satellite separations 
where the microlensing technique is optimized for the detections of moons. By setting an upper limit as the 
angle-average of the projected Hill radius and a lower limit as the half of the Einstein radius of the planet, we 
find that the microlensing method would be sensitive to moons with projected separations from the planet of 
0.05 AU < d p < 0.24 AU for a Jupiter-mass planet, 0.03 AU < d p < 0.17 AU for a Saturn-mass planet, and 
0.01 AU <d p < 0.08 AU for a Uranus-mass planet. We compare the characteristics of the moons to be detected 
by the microlensing and transit techniques. 
Subject headings: gravitational lensing 



1. INTRODUCTION 

All planets in our solar system except Mercury and Venus 
have moons. With the increasing number of discovered ex- 
trasolar planets, the existence of moons and their character- 
istics in these exoplanets emerge as new questions. Sev- 
eral methods to answe r these questions have been proposed. 
ISartoretti & Schneiderl {1999) pointed out that high-precision 
photometry of stars during planet transit can be used to detect 
extrasolar moons either by direct sate llite transit or per turba- 
tion in the timing of the planet transit. iBrown et al.l (1200 ll) ap- 
plied this method to the transit planet HD 209458b and placed 
upper limits on moons orbiting the planet by using the transit 
light curve obtained from Hubble Space Telescope observa- 
tions. 

In addition to the transit method, microlensing technique 
can also be used for the detections o f extrasolar moo n s. Thi s 
possibility was first mentioned by iBennett & Rhid (2002). 
They claimed that space-based lensing surveys with high pre- 
cision and cadence would be able to detect not only planets 
but also moons orbiting the plane ts. From the investigation of 
satellite-induced lensing signals, lHan & Hani d2002l) pointed 
out that detections of Earth-Moon like systems would be dif- 
ficult because the satellite signal would be seriously smeared 
out by severe finite source effect. However, moons with large 
masses may exist. From detailed investigation of the long- 
term dynamical stability of moons, Barn es & O'Brienl (2002) 
pointed out that Earth-like moons of Jovian planets could have 
stable orbits for long time scales. If such massive moons are 
common, it will be possible to detect them by using the mi- 
crolensing technique. 

In this paper, we investigate the characteristics of the lens- 
ing signals of Earth-like moons orbiting ice-giant planets. We 
investigate the variation of satellite signals depending on the 
locations of satellites and masses and locations of host plan- 
ets. We also determine the range of satellite separations where 
the microlensing technique is optimized for the detections of 



2. BASICS OF LENSING 

For the description of the lensing behavior produced by 
satellite systems, it is required to include at least three lens 
components of the host star, planet, and moon. For a multiple- 
lens system, the image mapping from the lens plane to the 
source plane is expressed as 



m k /M 
z-Zlm ' 



(1) 



where N is the number of the lens components, £ = £ + ir\, 
ZL.k = xlm + iyL.k, and z = x + iy are the complex notations of 
the source, lens, and image positions, respectively, z denotes 
the complex conjugate of z, are the masses of the indi- 
vidual lens components, M = £\ is the total mass of the 
system, and m^jM represent the mass fractions of the individ- 
ual lens components. Here all lengths are expressed in units 
of the Einstein radius that is related to the lens mass and the 
distances to the lens (Dl) and source (Ds) by 
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(2) 

Due to lensing, the image of the source star is split into mul- 
tiple fragments and the individual images are distorted. The 
fragmentation and distortion of the source image cause varia- 
tion of the source brightness. The lensing process conserves 
the source surface brightness, and thus the magnification of 
each image corresponds to the ratio between the areas of the 
image and source. For an infinitesimally small source, the 
magnification of each image is obtained by the Jacobian of 
the mapping equation, i.e. 



dz dz 
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For Galactic lensing events, the typical separations between 
images are of the order of 0.1 mas and thus the individual 
images cannot be resolved. However, events can be noticed 
by the variation of the source star flux where the total mag- 
nification corresponds to the sum of the magnifications of the 
individual images, i.e. A = 

One important difficulty in describing the lensing behav- 
ior of a multiple lens system is that the mapping equation 
is expressed in terms of the source position as a function of 
the image and lens positions. This implies that finding im- 
age positions for a given source position requires inversion 
of the mapping equation but the inversion is algebraically 
impossible for a multiple lens system. One way to obtain 
the image positions is expressing the mapping equation as 
a poly nomial in z and th en numerically solving the polyno- 
mial (I Witt & Maolll995l) . The advantage of this method is 
that it enables semi-analytic description of the lensing be- 
havior and saves computation ti me. Howev er, the order of 
polynomial increases as N 2 + 1 (iRhie 1 1 1991) and thus solv- 
ing the polynomial becomes difficult as the number of lens 
components increases. In this case, one can still obtain the 
magni fication patterns by using th e inverse ray-shooting tech- 
nique (ISchneider & Weissl Il986t iKavser. Refsdal & Stabelll 
Il986t IWambsganss. Paczvhski & Schneiderl 11990b . In this 
method, a large number of light rays are uniformly shot from 
the observer plane through the lens plane and then collected 
(binned) in the source plane. Then, the magnification pat- 
tern is obtained by the ratio of the surface brightness (i.e., the 
number of rays per unit area) on the source plane to that on 
the observer plane. Once the magnification pattern is con- 
structed, the light curve resulting from a particular source tra- 
jectory corresponds to the one-dimensional cut through the 
constructed magnification pattern. Although this method re- 
quires a large amount of computation time for the construction 
of detailed magnification patterns, it has an important advan- 
tage that the lensing behavior can be investigated regardless of 
the number of lenses. In addition, one can easily incorporate 
the finite source effect, which is important for the description 
of the perturbati ons caused by low- mass objects such as plan- 
ets and moons (Ben nett & RhidfT996h . Due to this reason, we 
use the ray-shooting method for the investigation of magnifi- 
cation patterns. 

Due to the small mass ratio of the planet and even smaller 
mass of the moon, the lensing light curve of an event pro- 
duced by a star having a planet with moons is well described 
by the single-lens light curve produced by the host star for 
most of the event duration. A short-duration perturbation oc- 
curs when the planet happens to be at the location of one of 
the two images of the source star produced by the host star 
JGaudT & Gould Tl997h . Since the moon is close to the planet, 
the moon can also perturb the image and produce an addi- 
tional anomaly. The position of the image-perturbing planet 
in the lens plane corresponds to the position of caustic in the 
source plane. In other words, perturbations occur when the 
source is located close to the caustic. The caustic represents 
the set of positions in the source plane at which the magnifi- 
cation of a point source event is infinite. For the binary lens 
case composed of the star and planet, there exist two sets of 
caustics. One is located very close to the star (central caustic) 
and the other is located away from the star (planetary caus- 
tic). Among the two perturbation regions around the individ- 
ual caustics, noticeable perturbations induced by the moon are 
expected only in the region around the planetary caustic. Two 
factors cause difficulties in finding satellite signatures around 



the central caustic region. First, the central perturbation pro- 
duced by the moon occurs in a very tiny region. As a result, 
the signature of the moon would be significantly washed out 
by the finite source effect. Second, the perturbation regions of 
the planet and the moon nearly coincide. Then, the anomaly 
in the lensing light curve would be dominated by that of the 
planet due to the overwhelming mass of the planet compared 
to the mass of the moon, making it even more difficult to iden- 
tify the satellite signature. We therefore focus on the perturba- 
tion region around the planetary caustic throughout the paper. 

The location of the planetary caustic is related to the star- 
planet separation by 

s c =s p (l-i) 2 , (4) 

where s p represents the position vector of the planet from the 
star and its length is normalized by the Einstein radius of the 
star. Then, the caustic is located on the planet side when the 
planet is outside the Einstein ring (s p > 1 .0), while it is lo- 
cated on the opposite side when the planet is inside the ring 
(s p < 1.0). The number of caustics also depends on the plan- 
etary separation and it is one when s p > 1.0 and two when 
s p < 1.0. The caustic is within the Einstein ring when the 
planetary separation is within the range of 0.6 < s p < 1.6. 
The caustic size, which is proportional to the chance of plan- 
etary perturbation, is maximized when the planet is in this 
region and thus this r egion is often called as 'lensing zone' 
dGould & Loeblfl992l) . As the separation departs from the 
Einstein radius, the caustic becomes smaller as oc s~ 2 for 

s p ^> 1.0 and oc s 2 for s p <C 1.0. In addition, the caustic size 
becomes smaller with the decrease of the planet/star mass ra- 
tio as oc qp 2 dHanll2006h . When the perturbation is produced 
by a planet located outside of the Einstein ring, the sign of the 
resulting anomaly in the lensing light curve is positive, im- 
plying that the magnification during the perturbation is higher 
than the corresponding magnification of the single lens event. 
On the other hand, if the perturbation is produced by a planet 
located outside of the ring, the sign of the anomaly is negative 
(lHan & Chandl2003h . 

3. MAGNIFICATION PATTERN 

We investigate the microlensing signals of an Earth-like 
moon around ice-giant planets. For this investigation, we con- 
struct magnification patterns of lens systems with physical pa- 
rameters adopted from those of typical galactic microlensing 
events c urrently being detected toward the galactic bulge di- 
rection (ISumi et al.ll2003h IUdalskill2003h . We assume that the 
planet-hosting star is located at a distance of = 6 kpc from 
the observer and has a mass of = 0.3 Mq. We also assume 
that the source star is located at £>s = 8 kpc, that corresponds 
to the distance to the Galactic center. Then the physical Ein- 
stein radius corresponding to the lens mass and distance is 
r% = Dl6e = 1 .9 AU. For the star-planet and planet-moon sep- 
arations, we test various combinations keeping in mind that 
the planet-moon separation should have an upper limit. This 
upper limit is usually set by the Hill radius which approxi- 
mates the gravitational sphere of influence of the planet in the 
face of the perturbation from the host star. The Hill radius is 
related to the semi-major axis, a, of the planet and the masses 
of the star, M±, and planet, m p , by 

/ m \ 1/3 
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FIG. 1. — Magnification patterns of lens systems composed of a star, 
a planet, and a moon. The masses of the individual lens components are 
0.3 Mq for the star, 10 Me for the planet, and 1.0 Me for the moon. Each 
map is centered at the center of the planetary caustic produced by the planet. 
The moons have a common position angle of<j> = 60° with respect to the star- 
planet axis, where the planet is located on the left. See Figure[2]for the geom- 
etry of the lens system. The labels above and on the right side represent the 
projected star-planet and planet-moon separations, respectively. The value 
in the parenthesis .v„ represents the star-planet separation in units of the Ein- 
stein radius corresponding to the mass of the host star, while s s represents the 
planet-moon separation normalized by the Einstein radius corresponding to 
the mass of the planet. Grey-scale is drawn such that brighter tone represents 
higher magnification. The panels blocked by thick solid lines represent the 
cases where the planet-moon separation is greater than the angle-average of 
the plant's Hill radius. The light curves resulting from the source trajectories 
marked by straight lines with arrows in the individual panels are presented in 
the corresponding panels of Fig. [3] 

Microlensing is only sensitive to the projected separation, 
while the 3 -dimensional separation is important for the orbital 
stability. We, therefore, set the angle-average of the projected 
Hill radius, i.e. y/2/3rn, as the upper limit of the planet-moon 
separation. Since the satellite signal is an additional perturba- 
tion to the planet-induced perturbation, moons would be de- 
tected for events where planets are detected. We, therefore, 
test planets located within the lensing zone of the host star. 
In physical units, this corresponds to 1.2 AU < dp < 3.0 AU, 
where d p is the projected star-planet separation. 

In Figure Q] we present the magnification patterns induced 
by a planet with moons of various projected separations from 
the planet. The planet has a mass of 10 Me, where Me is 
the mass of the Earth. Each map is centered at the center 
of the planetary caustic produced by the planet. The moons 
have a common position angle of <f> = 60° with respect to the 
star-planet axis, where the planet is located on the left. See 




FIG. 2. — Geometry of the lens system composed of a star, a planet, and 
a moon. The area in the box represents the region where the magnification 
pattern is presented in Figure[T] 

Figure [2] for the geometry of the lens system. The labels 
above and on the right side of the maps represent the pro- 
jected star-planet and planet-moon separations, respectively. 
The notations d s and s s represent the projected planet-moon 
separations expressed in physical units and in units of the Ein- 
stein radius corresponding to the mass of the planet, rE lP , re- 
spectively. Grey-scale is drawn such that brighter tone rep- 
resents higher magnification. The panels blocked by thick 
solid lines represent the cases where the planet-moon sepa- 
ration is greater than the angle-average of the projected Hill 
radius of the planet and thus moons are prohibited to reside. 
The light curves resulting from the source trajectories marked 
by straight lines with arrows in the individual panels are pre- 
sented in the corresponding panels of Figure [3] For the con- 
struction of light curves, we take the finite-source effect into 
consideration by assuming that the source star has a radius of 

From the magnification patterns and light curves, we find 
that non-negligible satellite signals occur when the planet- 
moon separation is similar to or greater than the Einstein ra- 
dius of the planet, i.e. s s > 1.0. One thing to be noted is 
that the satellite signal does not diminish with the increase 
of the planet-moon separation beyond the Einstein radius of 
the planet. This contrasts to the planetary signal that van- 
ishes when the planet is located well beyond the Einstein ra- 
dius of the star. This is because although the lensing effect of 
the planet on the moon rapidly decreases with the increase of 
the planet-moon separation beyond rB, p , the effect of the star 
on the moon remains. When, the planet-moon separation is 
substantially larger than the Einstein radius of the planet, the 
moon-induced perturbation forms at a separate region from 
the planet-induced perturbation region. In this case, the satel- 
lite signal on the light curve appears as a separate anomaly 
and thus it would be easily noticed. When s s ~ 1 .0, the per- 
turbations induced by the planet and satellite interfere each 
other, resulting in complex magnification patterns. Then, al- 
though it would be still possible to notice the satellite signal, 
it would be sometimes difficult to unambiguously identify the 
satellite signal. When the separation is substantially smaller 
than the planetary Einstein radius, the planet and moon be- 
have as if they are a single component. In this case, it would 
be difficult to notice the satellite signal. 

Another interesting trend of the satellite signal is that it 
tends to have the same sign as that of the planetary signal. 
This trend occurs because in most cases of identifiable moons 
with separations from the planet of s p > 1 .0, the lensing effect 
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FIG. 3. — Light curves of lensing events produced by lens systems com- 
posed of a star, a planet, and a moon. The lens system geometry and source 
trajectories responsible for the individual events are presented in the corre- 
sponding panels of Fig. [T] In each panel, the thick and thin curves represent 
the light curves resulting from lens systems with and without the moon, re- 
spectively. 
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FIG. 4. — Magnification patterns of lens systems obtained by using the 
binary superposition approximation. Notations are same as in Fig.fT] 



of the host star on the moon is bigger than the effect of the 
planet. Then, the sign of the satellite perturbation is mostly 
determined by the star-moon separation. The star-moon sep- 
aration is similar to the star-planet separation, and thus the 
sign of the planet-induced perturbation is same as that of the 
planet-induced perturbation. 

We note that although the planet and moon often reveal 
themselves as separate signals, characterizing them from the 
independent analysis of the individual signals would be diffi- 
cult. For some cases of triple lensing where the effect of the 
second body on the third body is negligible, it is possible to 
approximate the lensing behavior of the triple-lens system as 
the superposition of the two binary lens pairs composed of 
the first and second bodies and the first and third bodies. An 
example is the multiple-planetary system, wher e the lensing 
effect of a plane t to another planet is negligible dBozzall 19991 : 
lHan et all 1200 lh . However, for the case of the star-planet- 
satellite system, the effect of the planet on the moon is usually 
not negligible and thus the approximation of binary superpo- 
sition cannot be used for the analysis of the satellite signal. 
This can be seen from the comparison of the magnification 
patterns obtained by the exact triple-lensing formalism in Fig- 
ure [Hand the patterns obtained by using the binary superposi- 
tion approximation in Figure @] As expected, it is found that 
the difference in the magnification patterns becomes larger as 
the planet-moon separation decreases. 

Then, what will be the range of the satellite separation 
where the microlensing technique is optimized for the de- 
tections of moons. The lower limit of this range is set by 



the Einstein radius of the planet because moons with separa- 
tions substantially smaller than r E p are hard to be detected. 
Since r E p depends on the planet's mass, planets with different 
masses have different lower limits. The upper limit is set by 
the Hill radius because moons cannot reside beyond r^. The 
Hill radius depends not only on the planet mass but also on the 
star-planet separation. As a result, even planets with similar 
masses have different upper limits depending on where they 
are located in the system. 

In Figure |5J we present the optimal range of satellite sepa- 
rations as a function of the star-planet separation for planets 
with different masses. We note that the labels of the star- 
planet separation on the bottom axis (s p ) and the planet-moon 
separation on the left axis (s s ) are expressed in units of the 
Einstein radii of the star and planet, respectively. The la- 
bels are expressed also in physical units on the top and right 
axes, respectively. The individual panels are for planets with 
masses of m p = 300 M E , 100 M E , and 10 M E , which roughly 
corresponds to the masses of Jupiter, Saturn, and Uranus, re- 
spectively. In each panel, the light-shaded area represents the 
region of detectable satellites. The dark-shaded area repre- 
sents the region where the planet-moon separation is larger 
than the angle-average of the projected Hill radius and thus 
moons are prohibited to reside. The hatched area represents 
the region where the planet-moon separation is smaller than 
half of the Einstein radius of the planet and thus the satellite 
signal is hard to be detected. Although the range varies de- 
pending on the planet's position in the stellar system, we find 
that the microlensing method would be sensitive to moons 
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FIG. 5. — The optimal ranges of satellite separations where microlensing 
technique is sensitive to the detections of moons. The labels of the star- 
planet separation on the bottom axis (.v p ) and the planet-moon separation on 
the left axis (s s ) are expressed in units of the Einstein radii of the star and 
planet, respectively. The labels on the top and right axes are expressed in 
physical units. In each panel, the light-shaded area represents the region of 
detectable satellites. The dark-shade area represents the region where the 
planet-moon separation is larger than the angle-average of the projected Hill 
radius of the planet and the hatched area represents the region where the 
separation is smaller than half of the planetary Einstein radius. 

with separations from the planet of 0.05 AU < d p < 0.24 AU 
for a Jupiter-mass planet, 0.03 AU < d v < 0.17 AU for a 
Saturn-mass planet, and 0.01 AU < d p < 0.08 AU for a 
Uranus-mass planet. 

4. COMPARISON TO TRANSIT METHOD 

Due to the uniqueness of the microlensing method in detect- 
ing planets and their moons, the characteristics of the moons 
to be detected by the microlensing method will be different 
from those to be discovered by the transit method. Below, we 
list some of these differences. 

First, while the transit method can be used to search for 
moons of nearby stars, the microlensing method are sensitive 
to moons of remote stars. To meet the precision of photometry 
that is required to detect moons of extrasolar planets, the tar- 
get stars of transit searches should be bright and thus they are 
confined to the solar neighborhood. By contrast, microlens- 
ing searches are sensitive to stars anywhere along the line-of- 
sight toward the Galactic bulge. Therefore, the microlensing 
method can provide a sample of extrasolar moons distributed 
throughout the galaxy. 

Second, while the transit method is most sensitive to moons 
of close-in planets, the microlensing method is sensitive to 
moons of planets in the region beyond the 'snow line'. The 
snow line is the point in the protoplanetary disk beyond which 
the temperature is less than the condensation temperature of 




0.2 0.4 0.6 

planet-moon separation (r H ) 

FIG. 6. — The semi-major axis distributions of the moons of Jupiter, Saturn, 
and Uranus. Note that the planet-moon separations are expressed in units 
of the Hill radii of the individual planets. The inset in each panel shows 
the distribution for close-in satellites to the planet. The arrow represents the 
upper limit of the planet-satellite separation for the detection of moons by 
using the transit method. 

water dLecar et a l. 2006). Enhanced surface density of solids 
helps the formation of cores of giant planets and thus giant 
planets are thought to form in the region immediately beyond 
the snow line. The giant planets in our solar system, which 
are located in this region, have numerous moons; 63 known 
moons for Jupiter, 59 for Saturn, and 27 for Uranus. On the 
contrary, there might be few moons in close-in planets due to 
the strong tidal effect of the host stars as suggested by the two 
innermost planets of Mercury and Venus in our solar system. 

Third, while the transit method is sensitive to moons located 
close to their host planets, moons detectable by the microlens- 
ing method will have wide separations from the planets. In 
order to produce additional dips in transit light curves, moons 
should be located very close to their planets with separations 
equivalent to or less than the diameter of the host star. On the 
other hand, the sensitivity of microlensing method extends up 
to the Hill radius. For the case of giant planets in our solar 
system, the numbers of moons with separations larger than 
0.1r H are 54 (85.7% of the all known moons), 38 (64.4%), 
and 5 (18.5%) for Jupiter, Saturn, and Uranus, respectively. 
See the semi-major axis distributions of the moons of Jupiter, 
Saturn, and Uranus in Figure [6] In these planets, there also 
exist close-in moons with separations from the planets less 
than the diameter of the sun; 7 (11.1%), 19 (32.2%), and 17 
(63.0%) for Jupiter, Saturn, and Uranus, respectively. How- 
ever, it would be difficult to detect them by using the transit 
method because the star-planet separation is large and thus the 
probability of planet transit is very low. 
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5. CONCLUSION 

We investigated the characteristic of microlensing signals 
of Earth-like moons orbiting ice-giant planets. For this, we 
constructed magnification patterns of lens systems with vari- 
ous star-planet and planet-moon separations. From this inves- 
tigation, we found that non-negligible satellite signals occur 
when the planet-moon separation is similar to or greater than 
the Einstein radius of the planet. We found that the satellite 
signal does not diminish with the increase of the planet-moon 
separation beyond the Einstein radius of the planet unlike the 
planetary signal which vanishes when the planet is located 
well beyond the Einstein radius of the star. We also found 
that the satellite signal tends to have the same sign as that 
of the planetary signal. These tendencies are caused by the 
lensing effect of the star on the moon in addition to the effect 
of the planet. We determined the range of satellite separa- 
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detections of moons. By setting an upper limit as the angle- 
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ing method would be sensitive to moons with projected sep- 
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